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A  Monte  Carlo  Simulation  of  Guidance  Accuracy  Evaluation 


Wu  Liren 

(Beijing  Aerospace  Automatic  Control  Institute) 


Abstract:  It  is  necessary  to  carry  out  a  Monte  Carlo 

simulation  of  guidance  accuracy  evaluation  in  order  to 
solve  the  problem  of  evaluating  guidance  accuracy 
because  of  the  few  strategic  missile  flight  tests  that 
have  been  carried  out.  A  Monte  Carlo  simulation  method 
of  strategic  missile  guidance  accuracy  evaluation,  as 
well  as  several  important  problems  involved  in 
engineering  application  and  specific  methods  to  solve 
these  problems  are  introduced. 

Key  Words:  Strategic  missile,  guidance  error,  Monte 
Carlo  method,  mathematical  model 


1 .  Introduction 

Strategic  missile  accuracy  is  a  vitally  important  tactical 
and  technical  index.  To  evaluate  the  guidance  accuracy,  the 
conventional  method  is  to  conduct  a  statistical  evaluation  of  the 
missile  falling  point  during  a  flight  test,  which  requires  a 
certain  number  of  flight  tests.  In  reality,  however,  the  number 
of  missile  flight  tests  is  extremely  limited  due  to  their  very 
high  costs  and  the  restrictions  imposed  by  various  conditions  at 
test  sites.  Therefore,  a  strategic  missile  flight  test 
essentially  bears  the  character  of  accuracy  verification.  To 
overcome  the  difficulties  confronted  by  the  conventional  method 
in  this  area,  a  Monte  Carlo  simulation  of  strategic  missile 
guidance  accuracy  evaluation  has  been  proposed  and  is  regarded  as 
a  significant  technical  approach  in  solving  this  problem. 
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The  Monte  Carlo  method  is  an  experimental  mathematical 
method  with  which  a  statistical  check  can  be  carried  out  by  using 
random  numbers,  and  the  characteristic  statistical  values 
calculated  (such  as  mean  value,  variance,  probability,  etc.)  are 
taken  as  a  numerical  solution  to  a  particular  problem  to  be 
solved.  More  specifically,  this  is  a  method  of  reconstructing 
the  physical  process  of  a  practical  problem  under  research  by 
using  a  mathematical  probability  model  and  calculations  of  the 
characteristic  statistical  values  of  that  process.  With  this 
method,  the  model  established  can  be  tested  many  times,  and  based 
on  this  procedure,  the  test  data  are  processed  statistically  so 
as  to  derive  the  characteristics  of  the  foregoing  process  as  the 
statistically  estimated  values  of  the  process  parameters.  In 
accordance  with  the  amount  of  scattering  of  these  parameters,  the 
degree  of  approximation  in  solving  the  problem  can  be  determined 
in  the  sense  of  probability. 

The  Monte  Carlo  method  itself  is  well  developed.  With  the 
development  of  computers  in  the  past  few  decades,  this  method  has 
been  applied  more  and  more  widely  in  various  scientific  and 
technological  engineering  areas,  and  has  solved  a  large  number  of 
practical  problems.  This  paper  mainly  discusses  the  Monte  Carlo 
simulation  of  strategic  missile  guidance  evaluation,  as  well  as 
some  major  problems  that  this  kind  of  simulation  encounters  in 
its  engineering  applications  and  their  specific  solutions. 

2.  Basic  Aspects  of  Monte  Carlo  Simulation  of  Guidance  Accuracy 
Evaluation 

Simulating  strategic  missile  guidance  accuracy  with  the 
Monte  Carlo  method  proves  to  be  an  extremely  complicated 
numerical  simulation  test  of  engineering  systems .  It  is  supposed 
to  be  accomplished  in  three  steps,  namely  building  a  mathematical 
model,  conducting  a  simulation  test  and  analyzing  the  simulated 
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test  data.  This  procedure  involves  eight  major  aspects  listed  as 
follows : 

a)  Analysis  and  research  on  various  sources  of  errors 
affecting  the  guidance  system  accuracy,  as  well  as  collection  and 
accumulation  of  the  number  of  errors; 

b)  Constructing  probability  models  for  guidance  tool  errors 
and  guidance  method  errors; 

c)  Constructing  a  mathematical  model  of  the  trajectory  which 
is  used  to  describe  the  actual  missile  flight; 

d)  Confirmation  of  the  correctness  of  the  simulation  model 
established; 

e)  Selection  of  a  proper  random  variable  sampling  method; 

f)  Design  of  the  simulation  test  which  includes  correct 
generation  and  statistical  check  of  pseudorandom  numbers,  design 
of  a  method  of  minimizing  the  simulation  test  variance,  and 
determination  of  simulation  test  times  required  for  a  given 
accuracy; 

g)  Analysis  and  processing  of  simulation  test  data,  and 

h)  Computing  software  design,  debugging  and  verification  on 
a  computer  which  is  necessary  for  implementing  the  Monte  Carlo 
simulation  test. 

3 .  Method  of  Constructing  the  Mathematical  Model 

A  key  issue  in  evaluating  guidance  accuracy  with  Monte  Carlo 
simulation  is  how  to  correctly  establish  the  mathematical  model. 
Construction  of  the  mathematical  model  involves  two  major 
aspects : 

a)  Constructing  a  set  of  coupled  differential  equations 
which  can  precisely  describe  the  actual  missile  flight  movement, 
i.e.  a  mathematical  model  of  the  trajectory; 

b)  Fully  taking  advantage  of  the  information  from  existing 
ground  tests  and  flight  tests  so  as  to  collect  and  accumulate 
various  statistical  data  concerning  error  levels  and  build 
probability  models  for  differing  error  levels. 
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Among  other  things ,  the  construction  of  probability  models 
is  the  most  difficult  and  most  important  task  in  establishing  the 
mathematical  model.  The  probability  models  constructed  are 
required  not  only  to  be  able  to  reflect  the  reality  of  the  object 
simulated,  but  also  to  be  simple  and  easy  to  realize  so  that  the 
solution  coincides  with  the  probability  distribution  or 
mathematical  expectation  of  the  model  constructed.  Since  a 
simulation  test  without  a  certain  confidence  level  makes  no 
sense,  it  is  necessary  to  confirm  the  correctness  of  the 
simulation  model  constructed  before  the  simulation  test  can  be 
started . 

Based  on  an  analysis  of  the  various  levels  of  errors  which 
affect  guidance  accuracy,  two  methods  of  building  error  level 
probability  models  are  introduced. 


3.1  Analysis  of  Various  Error  Levels  Affecting  Guidance  Accuracy 

A  ballistic  missile  system  may  have  two  kinds  of  error 
sources :  guidance  errors  and  non-guidance  errors .  The  inertial 
guidance  system  consists  of  an  attitude  reference,  an 
accelerometer,  a  computer  and  a  control  system.  The  guidance 
errors  can  be  further  subdivided  into  guidance  tool  errors  and 
guidance  method  errors. 

The  problems  which  affect  the  guidance  tool  errors  include: 
initial  alignment  errors  (leveling,  pointing)  of  a  stable 
platform  before  the  launch;  errors  in  missile  velocity 
measurement  with  the  accelerometer;  missile  attitude  measurement 
error  caused  by  gyrodrift;  in-missile  computer  errors,  etc. 

The  problems  which  affect  the  guidance  method  errors 
include:  missile  launch  quality  deviation;  deviation  in  engine 
propellant  agent  per-second  consumption;  engine  thrust  ratio 
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deviation;  engine  thrust  line  inclination;  deviation  caused  by- 
wind  interference  and  atmospheric  parameter  changes ,  etc . 

In  terms  of  kinematic  trajectory  equation  calculations,  the 
foregoing  problems  can  be  subdivided  into  three  categories  as 
follows : 

a)  There  is  no  exact  knowledge  about  the  number  of  constants 
in  the  kinematic  equation,  such  as  take-off  quality,  engine 
thrust  line  inclination; 

b)  The  actual  parameter  change  regularities  do  not  conform 
to  the  design,  such  as  engine  propellant  agent  per-second 
consumption  deviation,  variation  of  atmospheric  parameters  with 
altitude,  etc.; 

c)  Parameters  which  are  not  included  in  the  kinematic 
equation,  such  as  inertial  instrument  errors,  non-uniformity  of 
wind  and  atmospheric  parameters,  etc. 

Any  of  the  foregoing  problems  has  a  random  character,  i.e. 
their  effect  on  any  particular  in-flight  missile  cannot  be 
predicted  in  advance.  Every  error,  although  it  has  an  uncertain 
value,  abides  by  a  given  statistical  distribution  and  therefore 
its  effect  on  guidance  errors  can  be  calculated  as  long  as  its 
distribution  is  known. 

3.2  Method  of  Constructing  a  Probability  Model 

When  the  original  data  on  various  errors  which  affect  the 
guidance  accuracy  are  gathered  and  properly  processed,  a  set  of 
observation  data  xl ,  x2...xN  which  reflect  a  particular  random 
process  of  a  particular  error  can  be  acquired.  If  these 
observation  data  are  independent  of  each  other,  the  model  can  be 
constructed  by  selecting  the  random  variable  probability 
distribution,  or  otherwise  if  they  are  dependent,  correlated  to 
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each  other;  the  model  should  be  built  with  the  time  sequence 
analysis  method. 

3.2.1  Model  Construction  Based  on  Selection  of  Random  Variable 
Probability  Distribution 

As  mentioned  above,  when  the  observation  data  which  reflect 
a  random  variable  of  a  given  error  are  independent,  a  random 
variable  probability  distribution  can  be  selected  to  construct  a 
probability  model.  To  achieve  this  goal,  the  usual  practice  is 
to  select  a  proper  probability  distribution  type  before  carrying 
out  a  parameter  estimation. 

3 . 2 . 1 . 1  Method  of  Selecting  Probability  Distribution  Type 

Based  on  a  priori  knowledge  including  an  understanding  of 
the  particular  random  variable  and  data  possessed,  it  can  be 
decided  which  distribution  should  be  adopted  or  which 
distribution  should  be  rejected.  In  this  case,  data  play  a 
decisive  role  in  selecting  the  probability  distribution.  With 
adequate  data,  it  is  possible  to  select  a  probability 
distribution  which  conforms  with  the  actual  situation  of  the 
system,  and  also  to  check  the  correctness  of  the  probability 
distribution  selected. 

The  following  description  is  given  for  some  specific  methods 
which  can  be  used  in  selecting  an  appropriate  probability 
distribution  type  with  given  data: 

3. 2. 1.1.1  Histogram  Method 

A  histogram  is  a  graphic  description  of  data  evaluation 
distribution,  and  also  an  approximation  of  the  random  variable 
density  function.  By  estimating  the  shape  of  a  density  curve 
with  a  histogram,  it  can  be  roughly  seen  which  type  of 
probability  distribution  it  looks  like. 
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3. 2. 1.1. 2  Probability  Diagram  Method 


The  probability  diagram  method  is  used  to  estimate  the 
theoretical  distribution  function  directly  with  the  empirical 
distribution  function  FN(x) .  The  method  generally  used  in 
engineering  is  the  probability  paper  method.  With  the 
probability  paper  method,  some  theoretical  distributions  can  be 
put  in  a  given  coordinate  in  straight  lines,  and  in  the  meantime, 
an  empirical  distribution  function  curve  is  also  plotted  in  the 
same  coordinate  system  for  a  comparison  to  find  out  their 
resemblance  and  difference.  For  instance,  the  probability 
diagram  of  a  particular  datum  can  be  drawn  in  a  normal  state 
probability  paper,  and  if  the  datum  is  distributed  in  a  straight 
line,  it  is  considered  to  abide  by  the  normal  state  distribution. 
In  this  case,  many  kinds  of  probability  papers  can  be  prepared  in 
advance,  including  a  normal  state  probability  paper,  an 
algorithmic  normal  probability  paper,  a  Weber  probability  paper, 
etc . 

3 . 2 . 1 . 1 . 3  Point  Estimation  Method 

With  this  method,  sample  variance  S2  and  sample  mean  value 
can  be  calculated  in  accordance  with  a  given  sample  datum  and 
further,  the  estimated  value  ,7  of  the  variance  coefficient  5c 
can  be  calculated  on  the  basis  of  the  following  equation 


.x 


There  are  different  variance  coefficient  Sc  evaluation 
ranges  in  response  to  different  probability  distributions,  which 
can  be  found  in  a  table.  It  can  be  judged  then,  based  on  the 
value  calculated,  which  probability  distribution  is  obeyed  by  the 
foregoing  sample  datum. 
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3 . 2 . 1 . 2 .  Parameter  Estimation  Method 


After  a  probability  distribution  type  is  selected,  the 
parameters  of  that  distribution  should  be  estimated.  The 
commonly  used  methods  of  estimating  such  parameters  are  the 

moment  estimation  method  and  the  maximum  similarity  estimation 
method . 

3 . 2 . 1 . 2 . 1  The  Moment  Estimation  Method 

Since  Xi,  x2...xN  are  known  as  a  sample  of  the  random 
variable  X,  the  sample  k-order  original  point  moment  can  be  used 

as  the  estimation  of  the  overall  k-order  original  point  moment  as 
shown  in  the  following  equation 


2  ^or  instance,  it  is  known  that  the  random  variable  x~N(}i, 
a2),  the  moment  estimated  amounts  of  4  and  a2  are 


<7  =/i. 


3. 2. 1.2. 2  The  Maximum  Similarity  Estimation  Method 

The  basic  idea  of  the  maximum  similarity  estimation  method 
signifies  that  if  the  maximum  value  of  the  similarity  function 
ln(Xi,  x2...xn,  6i,  02...eM)  is  selected  at  0  ,  0 ,,  $ 

(i.e.  the  probability  of  their  generation  is  the  maximum 

parameters  s  ?,  -  \ 

0r  0,,  o i*  then  ^  ^  ...  q  can  be 


can  be 


8 


taken  as  the  maximum  similarity  estimated  value  of  the  overall 
parameter  0W  92,...,eM. 

As  for  the  normal  state  distribution,  its  maximum  estimated 
similarity  amount  is  identical  to  the  moment  estimated  amount, 
while  other  distributions  can  be  directly  calculated  by  using 
many  maximum  similarity  estimation  formulas  for  common 
probability  distribution  parameters,  which  are  listed  in 
statistics  books,  as  is  done  in  practical  computations. 

3.2.2  Model  Construction  Based  on  Time  Series  Analysis  Method 

3. 2. 2.1.  The  Time  Series  Analysis  Method 

During  quantitative  measurement  of  one  variable  or  a  group 
of  variables  X(t)  derived  from  a  production  or  research  process, 
a  discrete  digitalized  ordered  number  ensemble  x(ti )  ,x( t2) .  .  .x( tH) 
with  time  t  as  an  independent  variable  can  be  obtained  at  the 
moment  tx<t2< .  .  . <tN,  which  is  generally  referred  to  as  a  digital 
time  series .  In  practical  problems ,  the  independent  variable  t 
represents  not  only  time,  but  also  various  physical  implications, 
including  length,  temperature,  velocity  or  other  physical 
quantities  of  monotonically  increasing  evaluation. 

Generally  speaking,  it  is  extremely  difficult  to  express  a 
time  series  with  a  completely  certain  function  or  function  group. 
In  most  cases,  its  statistical  law  is  displayed,  and  therefore 
its  evaluation  law  can  be  described  statistically  through  a 
probability  distribution  function  or  function  group.  This  kind 
of  time  series  is  referred  to  as  a  random  time  series,  or  a 
random  process,  as  it  is  called  in  the  probability  theory.  Most 
time  series  can  be  regarded  as  a  reality,  or  a  sample  function  of 
the  random  process. 
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If  time  series  data  are  not  independent  but  correlated  to 
each  other,  this  correlation  must  be  taken  into  account  in  model 
construction,  i.e.  model  construction  should  be  based  on  the  time 
series  analysis  method  instead  of  the  probability  distribution 
selection  method. 

The  time  series  analysis  method  is  a  statistical 
mathematical  method,  which  is  used  for  time  series  analysis  and 
research.  The  major  time  series  analysis  methods  are  all  based 
on  the  assumption  that  data  samples  are  all  derived  from  a  smooth 
random  process,  i.e.  their  statistical  characteristics  such  as 
mean  values,  variances  and  correlation  functions  do  not  change 
with  the  passage  of  the  original  time  point.  In  reality, 
however,  most  random  processes  are  not  smooth,  and  therefore,  the 
non-smooth  time  series  must  initially  be  smoothed. 

The  parameter  method  and  the  difference  method  are  commonly 
used  as  smoothing  methods  for  non-smooth  time  series .  In  this 
paper,  our  discussion  is  devoted  to  the  parameter  method  alone. 

The  non-smooth  time  series  x(t)  is  usually  described  with 
the  following  addition  model  (other  non-addition  models  can  be 
properly  converted  to  addition  models ) : 

x(t)=d(t)+  n(t)  (4) 

where  tj(t)  is  assumed  to  be  a  time  series;  d(t)  is  a  function 
term  indicating  the  time  series  x(t)  with  an  apparent  tendency, 
which  is  presumably  formed  through  super imposition  of  the 
following  two  functions: 

d(t)=f(t)+p(t)  (5) 

where  f(t)  is  a  main  value  function  term  indicating  a  long-term 
change  tendency  of  the  time  series;  p(t)  is  a  periodic  function 
term  which  represents  some  periodic  change  factors  affecting  the 
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time  series.  Their  identification  and  extraction  methods  [3] 
briefly  introduced  in  the  following  section. 


are 


3. 2. 2. 1.1  Method  of  Identifying  and  Extracting  Main  Value 
Function  Term  f(t) 


The  following  functions 

J(t)  -a0  +*  /  +  *,/'  +a  , 


i*  +  V  +■*,/ 


+  Of  /  -f  Of  / -  -f  *  f  2  ,  ,  ■ 

6  7  V  +a|0ln/ 

are  used  as  an  approximation  of  the  actual  main  value  function 
term  so  as  to  fit  into  the  given  time  series  x(t).  In  other 
words,  in  the  the  method  of  least  squares,  a  computer  is  used  to 
identify  and  select  the  function  f(t),  estimate  the  regression 
coefficient  a±  (i=0,  1,  2... 10),  and  finally  give  the  form  of  the 
main  value  function  term  f(t)  by  using  the  gradual  regression 
analysis  method  and  under  the  checked  values  Fx  and  F2  of  the 
given  round-off  forecast  factor  F.  If  the  estimated  regression 
coefficient  a±  is  zero  in  any  case,  it  can  be  accepted  then  that 
there  is  no  main  value  function  term  f(t). 

3. 2. 2. 1.2  Method  of  Identifying  and  Extracting  Periodic  Function 
Term  p ( t ) 

When  the  main  value  function  term  f(t)  is  removed  from  the 
given  time  series  x(t),  an  ordered  number  ensemble  is  obtained  as 


Xi(t)=X(t)-f (t) 


Consider  the  periodic  series  model 


*,(0  =a0  +  Z  (a.-coseo./  + /J.sinco./)  +  rj(t) 

i «  1 

where  a0,  air  $±r  6^  i=l,2,...,l)  and  1  are  all  parameters  to  be 
identified  and  estimated. 


Now,  the  periodic  function  term  p(t)  can  be  identified  and 
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extracted  with  the  method  of  periodic  diagram  analysis  in  the 
least  square  method. 


After  the  calculated  function  term  d(t)  is  removed  from  a 
given  non-smooth  time  series  x(t),  the  time  series  ti(t)  can  be 
derived  as 

n(t)=x(t)-d(t)=x(t)-f (t)-p(t)  (6) 


Obviously,  ti(t)  is  a  smooth  time  series. 

3 . 2 . 2 . 2  Construction  of  Smooth  Linear  Model 

As  soon  as  the  smooth  time  series  n(t)  is  derived,  it  is 
necessary  to  determine  what  kind  of  smooth  linear  model  should  be 
adopted  through  model  identification  method.  The  commonly  used 
smooth  linear  models  are  the  self-regression  model  AR(p),  sliding 
average  model  MA(q)  and  self-regression-sliding  average  model 
ARMA(p,  q). 

The  time  series  model  is  established  in  the  order  of  model 
identification,  its  type  and  order  number  determination  and  its 
parameter  estimation. 

The  AR(p)  model  parameters  are  estimated  with  a  linear 
estimation  method.  In  this  case,  the  least  square  method  can  be 
used  to  derive  the  estimated  model  parameters  derived  from  a 
solution  of  the  standard  equation.  As  for  the  parameters  of  the 
MA(q)  and  ARMA(p,  q)  models,  their  estimation  appears  much  more 
complex  and  requires  the  application  of  the  non-linear  least 
square  method . 

It  is  to  be  noted  that  the  foregoing  time  series  analysis 
methods  of  model  construction  are  applicable  only  in  the  case 
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when  the  data  concerned  obey  the  normal  state  of  distribution, 
while  the  input  process  of  an  actual  system  is  often  a  non-normal 
state  process.  Thereby,  data  should  be  converted  from  a  non¬ 
normal  distribution  to  a  normal  distribution  before  the  foregoing 
methods  can  be  used.  Specifically,  the  conversion  methods 
include  the  algorithmic  conversion  method  and  the  normal 
conversion  method  based  on  the  empirical  distribution  function 
[2]. 

4.  Random  Number  Generation  and  Verification  Methods 

The  foundation  of  Monte  Carlo  simulation  is  that  a  digital 
computer  can  rapidly  generate  random  numbers  with  ideal 
statistical  properties,  which  are  evenly  distributed  in  the  [0, 

1]  interval.  With  random  numbers,  the  random  number  series  with 
different  distribution  regularities  required  from  the  simulation 
can  be  derived  through  conversion,  sampling,  etc.  The  quality  of 
random  series  has  a  direct  effect  on  the  confidence  level  of  the 
simulation  result. 

Generally,  random  numbers  are  generated  by  using 
mathematical  methods,  such  as  the  multiplication  congruence 
method,  addition  congruence  method,  mixed  congruence  method  and 
quadratic  method.  Procedures  are  readily  available  for  all  these 
methods  to  be  selected. 

As  for  the  pseudorandom  series  generated  on  a  computer  with 
mathematical  methods,  they  have  no  real  random  character.  These 
pseudorandom  numbers  must  be  checked  statistically  before  being 
used  as  required  random  numbers,  so  that  there  will  not  be  any 
large  systematic  error  in  the  Monte  Carlo  simulation  result. 

There  are  many  random  number  statistical  verification 
methods,  which  include  parameter  checking,  uniformity  checking, 
independence  checking,  continuity  checking  and  combination 
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regular  checking,  etc.  It  is  necessary  to  select,  based  on  the 
specific  problem  simulated,  an  appropriate  verification  method, 
and  determine  its  proper  marked  level  and  rejection  domain. 

5.  Falling  Point  Accuracy  Evaluation  Method 

Suppose  n  Monte  Carlo  simulation  tests  have  been  conducted, 
and  the  observation  values  of  falling  point  deviation  caused  by  a 
guidance  error  have  been  calculated  as 


then  the  falling  point  accuracy  and  accuracy  evaluation  can  be 
carried  out  over  the  values  of  that  sample  group  as  described  in 
the  following  section. 

5.1  Sample  Characteristic  Value  Calculations 
Sample  mean  value 

Sample  variance 
Sample  standard  difference 


n  .  ,  • 


s  =  'J~sr 


(?) 


5.2  Estimation  of  Confidence  Interval  between  Mean  Value  and 
Variance  With  a  Given  Confidence  Level 


Generally,  missile  falling  point  deviation  obeys  the  normal 
state  distribution  (n,  ct2).  The  following  discussion  contributes 
to  the  estimation  of  the  confidence  interval  of  \i,  a1. 

5.2.1  Estimation  of  Mean  Value  Confidence  Interval 

5. 2. 1.1  With  a  Known  Variance  (<7:  =  ,r02) 
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Under  a  given  confidence  level  (1-a),  the  confidence 
interval  of  mean  value  n  is: 


Ly  —  u 


x  +  u  • 


(8) 


I -  »  t  u  -  -7= 

V/7  *  v7T/ 

where  n — sample  capacity; 

s'  sample  mean  value; 

u0— a— horizontal  bilateral  fractile  of  standard  normal 
state  distribution,  which  is  listed  in  the  related  table. 


5 . 2 . 1 . 2  With  an  Unknown  Variance 


With  a  given  confidence  level  (1-a),  the  confidence  interval 
of  mean  value  is: 


-  i 


x  +  t 


(9) 


where  t*  is  a — horizontal  bilateral  azimuth  number  of  t 
distribution  with  degree  of  freedom  (n-1). 

5.2.2  Estimation  of  Variance  Confidence  Interval 


5. 2. 2.1  With  a  Known  Mean  Value  fi,  (h=ji0) 

With  a  given  confidence  level  (1-a),  the  confidence  interval  of 
variance  a2  is : 


X 


X(*  -n0y 


i 


(10) 


where  x n  is  q — horizontal  upper  side  fractile  of  x2  distribution 
with  a  degree  of  freedom  n;  q  respectively  represents  a/2  and 
1-a/ 2 . 


5 . 2 . 2 . 2  With  an  Unknown  Mean  Value  pi 
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With  a  given  confidence  level  (1-a),  the  confidence  interval 
of  variance  o2  is: 


~  1  )S2 


(OjzJJjS 


X* 


7*  n  ~~  1 


'  I 


;) 


(11) 


6 .  Method  of  Determining  Simulation  Test  Times  Required  from  a 
Given  Accuracy 

Since  Monte  Carlo  simulation  for  strategic  missile  guidance 
accuracy  evaluation  requires  a  host  of  tests,  and  the  trajectory 
calculations  are  time-consuming,  a  heavy  volume  of  computer  time 
has  to  be  spent  to  ensure  high  precision  simulation  calculations. 
Under  such  a  scenario,  research  should  be  carried  out  on  how  to 
determine  the  minimum  simulation  test  times  nA  without 
sacrificing  the  required  simulation  accuracy  in  solving  specific 
engineering  problems . 

The  following  equation  reflects  the  relationships  among  the 
confidence  level  (1-a)  of  the  estimated  value,  the  confidence 
interval  2  calculated  and  the  simulation  times  n  realized. 


P{\x  —  -x|  <e}  -  1  —  a 


(12) 


A  formula  is  derived  from  Eq.  (12)  which  indicates  that  the 
test  times  nA  needed  for  ensuring  a  given  accuracy  rA  in 
simulation  calculations  can  be  derived  as  follows: 


n  A  —  min 


.  S 

•  ) -  t 

v  Yl  2 


(13) 


where  is  the  half  length  of  mean  value  n  confidence 

interval  under  an  unknown  variance  a2,  i.e.  £  ; 

rA  is  an  absolute  accuracy  of  the  confidence  interval. 
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which  is  defined  as  the  half  length  of  the  given  confidence 
interval;  , 

t„  is  a — horizontal  bilateral  fractile  of  t  distribution 
with  degree  of  freedom  (n-1),  which  can  be  found  in  the  related 
table. 

7 .  Method  of  Minimizing  Test  Variance 

7.1  Significance  of  Minimizing  Test  Variance 

When  the  Monte  Carlo  method  is  used  to  carry  out  simulation 
calculations,  its  method  error  e  can  be  written  as 


where  r-  is  the  given  simulation  error; 

c  is  the  standard  difference  of  the  random  variable 
selected; 

n  is  test  times. 

Equation  (14)  shows  that  there  are  two  approaches  to 
increase  simulation  accuracy  with  the  given  confidence  level 
(1-a):  one  is  to  increase  test  times  n,  while  the  other  is  to 
minimize  variance  ct2. 

1 

As  £  is  directly  proportional  to  ,  a  10-fold  increase 

of  simulation  accuracy  requires  a  100-fold  increase  in  test 
times,  which  inevitably  results  in  a  considerable  increase  in 
computer  time.  On  the  other  hand,  since  s  is  in  direct 
proportion  to  a,  minimizing  variance  a2  will  bring  about  a 
striking  positive  result.  Theory  and  practice  suggest  that  a 
combination  of  various  variance  minimization  methods  and  skills 
can  raise  the  simulation  accuracy  by  10  or  even  100  times. 
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7.2  Methods  of  Minimizing  Test  Variance 


There  are  many  methods  which  can  be  used  to  minimize  test 
variance.  Additionally,  different  methods  can  be  used  for  the 
same  given  problem.  Here  the  skill  of  minimizing  test  variance, 
the  dual  variable  method,  is  introduced;  it  is  well-suited  to 
engineering  problems. 


The  major  idea  of  the  dual  variable  method  is  as  follows: 
Suppose  xf  and  yf  are  determined  from  distribution  f(x,  y) 
through  sampling,  with  a  mathematical  expectation  of  0(x£,  yf) . 
x'f  and  y'f  are  also  determined  from  the  same  distribution  f(x, 
y)  through  sampling,  whose  mathematical  expectation  is  0(x£,y'f). 
If  0(xf,  yf)  and  0(x'£,  y'£)  are  in  negative  correlation,  this 
pair  of  estimations  with  reverse  variation  but  mutual 
compensation  are  referred  to  as  a  dual  variable.  Therefore, 


can  be  taken  as  the  estimated  value  of  0,  which  is  a  non-deviated 
estimation.  The  variance  of  the  estimated  value  $  is 

y ) |/ (-Y ,  >•)} 

When  the  correlation  coefficient  p=-l,  the  variance  of  the 
estimated  value  |  as  a  dual  variable  sampling  skill  is  equal  to 
zero . 


Thus,  the  dual  variable  method  is  virtually:  to  find  one  (or 
several)  random  variable(s)  with  the  same  mathematical 
expectation  value  as  that  of  the  original  random  variable  (the 
mathematical  expectation  value  of  this  random  variable  is  the 
quantity  to  be  calculated),  and  try  as  much  as  possible  to  place 
them  in  negative  correlation.  Based  on  this,  the  linear 
combination  of  the  two  (or  more)  random  variables  is  taken  as  the 
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4 


A 


estimated  value.  By  using  this  method,  the  test  variance  can  be 
minimized  remarkably. 

This  paper  was  received  on  March  18,  1995. 
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